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Abstract 
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NMR is emerging as a valuable testbed for the investigation of foundational ques- 
tions in quantum mechanics. The present paper outlines the preparation of a class 
of mixed states, called pseudo-pure states, that emulate pure quantum states in the 
highly mixed environment typically used to describe solution-state NMR samples. It 
also describes the NMR observation of spinor behavior in spin 1/2 nuclei, the simula- 
tion of wave function collapse using a magnetic field gradient, the creation of entangled 
(or Bell) pseudo-pure states, and a brief discussion of quantum computing logic gates, 
including the Quantum Fourier Transform. These experiments show that liquid-state 
NMR can be used to demonstrate quantum dynamics at a level suitable for laboratory 
exercises. 
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1 Introduction 

The fundamental physics of NMR is again, 50 years after its discovery the subject of much 
discussion. The impetus behind this recent interest is the dramatic potential of quantum 
information processing (QIP) [9991], particularly quantum computing, along with the real- 
ization that liquid-state NMR provides an experimentally accessible testbed for developing 
and demonstrating these new ideas [02-91, 12, 1999, 1]. 

Most descriptions of quantum information processors have focused on the preparation, 
manipulation, and measurement of a single quantum system in a pure state. The appli- 
cability of NMR to QIP is somewhat surprising because, at finite temperatures, the spins 
constitute a highly mixed state, as opposed to the preferred pure state. However, NMR 
technology applied to the mixed state ensemble of spins (the liquid sample) does offer sev- 
eral advantages. Decoherence, which plays a detrimental role in the storage of quantum 
information, is conveniently long (on the order of seconds) in a typical solution sample, and 
it acts on the system by attenuating the elements of the density matrix and rarely mixes 
them. NMR spectrometers allow for precise control of the spin system via the application 
of arbitrary sequences of RF excitations, permitting the implementation of unitary transfor- 
mations on the spins. Effective non-unitary transformations are also possible using magnetic 
field gradients. The gradient produces a spatially varying phase throughout the sample, and 
since the detection over the sample is essentially a sum over all the spins, phase cancellations 
from spins in distinct positions occur. These characteristics of NMR enable the creation of 
a class of mixed states, called pseudo-pure states, which transform identically to a quantum 
system in a pure state [9993]. 

NMR does have several noteworthy disadvantages. A single density matrix cannot be 
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associated with a unique microscopic picture of the sample, and the close proximity of the 
spins prevents the study of non-local effects. Additionally the preparation of pseudo-pure 
states from the high temperature equilibrium state in solution NMR entails an exponential 
loss in polarization. [39991] 

In this paper, we review the results of a number of simple NMR experiments demon- 
strating interesting quantum dynamics. The experiments illustrate spinor behavior under 
rotations, the creation and validation of pseudo-pure states, their transformation into "en- 
tangled" states, and the simulation of wave function collapse via gradients. Additionally, 
the implementations of basic quantum logic gates are described, along with the Quantum 
Fourier Transform. 
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2 The Spin System 



The experiments were performed on the two-spin hetero nuclear spin system, 13 C-labeled 
chloroform ( 13 CHC1 3 ), thereby eliminating the use of shaped RF pulses. The 13 C (I) and 
the 1 H (S) nuclei interact via weak scalar coupling, and the Hamiltonian for this system is 
written as 



where uji and u>s are the Larmor frequencies of the 13 C and 1 H spins respectively and J « 
\ujj — u)s\ is the scalar coupling constant. 

In the standard model of quantum computation, the quantum system is described by a 
pure state. However, liquid-state NMR samples at room temperature are in highly mixed 
states, requiring the state of the system to be described by the density operator. In a liquid 
sample, the inter-molecular interactions are, for most practical purposes, averaged to zero 
so that only interactions within a molecule are observable; in other words, the sample can 
be thought of as an ensemble of quantum processors, each permitting quantum coherence 
within but not between molecules. For the purposes of this paper, the large density matrix 
of size 2 N x 2 N , where N is the number of spins in the sample, may be replaced by a much 
smaller density matrix of size 2 n x 2 n , where n is the number of distinguishable spin-| nuclei 
in the molecule. In the high temperature regime (e = ~ O(l0~ 6 )) the equilibrium 

density operator for the ensemble is 



H = wil z + uj s S z + 2irJI z S : 



(1) 





where the relative value of the gyromagnetic ratios is js/li ~ 4. 
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From the above, it is clear that at room temperature a spin system cannot be prepared 
in a pure state. However, it is possible to prepare a pseudo-pure state that transforms like 
a pure state. Also, notice that since the identity part of the density operator is invariant 
under unitary transformations, it is the deviation part of the density operator, that holds the 
information on the spin dynamics. Henceforth in this paper, the deviation density matrix 
will be simply referred to as the density matrix. The density operator is often written in the 
product operator basis formed by the direct product of individual spin operators[451999, 1]. 
The product operator technique is used throughout this paper to express the dynamics of 
the spin system. Furthermore, if n spins are coupled to one another, any arbitrary unitary 
operation can be composed from a series of RF pulses, chemical shift evolution and scalar 
coupling evolutions. [64999, 1] 
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3 Preparation of Pseudo-Pure States 

Before describing the creation of the pseudo-pure state, it is convenient to begin with a 
system of equal spin populations. This is achieved by applying the pulse sequence 



'71' 
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■■ s (±) 


'71' 


I,S 


.2. 


Uj) 


.2. 


y Uj) 
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[grad(z)] 



(3) 



to the equilibrium density matrix, resulting in 



(4) 



which has a balanced spin population. Because the eigenvalue structure of this density matrix 
is different from that of thermal equilibrium, there is no unitary transformation which could 
transform one to the other. The non-unitary gradient (where the non-unitarity refers to 
the spatial average over the phases created by the gradient) at the end of the above pulse 
sequence makes this transformation possible. Figure |l] shows a spectrum obtained after 
applying this sequence. 

Since the identity part of the equalized density matrix is unaffected by unitary transfor- 
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mations and undetectable in NMR, only the deviation density matrix, 
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which represents the excess magnetization aligned with the external magnetic field, is of 
interest. The above matrix representation has been made in the eigenbasis of the unperturbed 
Hamiltonian, and here the rows and columns have been labeled explicitly to avoid ambiguity. 
In the subsequent matrix expressions, the labels will be dropped. 

QIP requires the ability to create and manipulate pure states. NMR systems, however, 
are in a highly mixed state at thermal equilibrium. While single spin manipulation is not 
feasible in NMR, Cory et. al. [2, 3999, 5] have developed a technique by which the equilibrium 
state is turned into a pseudo-pure state. Such a state can be shown to transform identically 
to a true pure state as follows: according to the rules of quantum mechanics, a unitary 
transformation U maps the density matrix p to ft = UplA. Thus an iV-spin density matrix 
of the form p = (1 + (ip\) / 2 N is mapped to 



i + 



2 N 



(6) 



This shows that the underlying spinor \ip) is transformed one-sidedly by U just as a spinor 
which describes a pure state would be. 
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After equalizing the spin population from the thermal equilibrium state (eq. (5)), the 
application of 



"7T" 


i,s 




'n~ 


A. 


X 




.6. 



i,s 



[grad(z)} 



(7) 



results in the pseudo-pure state (neglecting the initial identity component) 



'll + J^(I z + S z + 2I z S z 



( \ 
10 









Figure || shows a series of spectra confirming the preparation of a pseudo-pure state. 
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4 Spinor Behavior 

Particles of half-integral spin have the curious property that when rotated by 2n, their wave 
functions change sign while a An rotation returns their phase factors to their original value. 
The change in the sign of the wavefunction is not observable for a single particle, but it can 
be seen through an interference effect with a second "reference spin." Spinor behavior, as 
this effect is called, was first experimentally measured using neutron interferometry [781999, 
1] and later using NMR interferometry [99991]. 

The following simple experiment describes how the spinor behavior can be seen in chlo- 
roform, where the spinor behavior of 13 C is correlated with the 1 H nuclei as a multiplicative 
phase factor. Consider the unitary transformation 



U 



10 

cos(f) -sin(f 

10 

sin(f) cos(f) 



\ 



) 



(9) 



As explained in section 6, this can be viewed as a rotation by <fi of the 13 C conditional on 
the 1 H being in the down state. This can be implemented via the pulse sequence 
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— X 



(10) 



Application of this pulse sequence to the state 2I Z S X , where the spinor behavior of the I-spin 
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is revealed by its correlation to the S-spin, results in 



2 coa{<f>/2)I t S x + 2 sm(4>/2)I x S x . 



It can be clearly seen that when <p = 2n the initial state gains a minus sign, but when 
= 47r the state returns to its initial value. The state 2I Z S X is made observable under 
the evolution of the internal hamiltonian previously defined and can be created from the 
equalized equilibrium state (eq. |j) using the sequence 



2 



[grad(z)\ 



2 



1 

2J 



(12) 



Figure H shows the spectra for several values of = 0, 2n, and 47r. 
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5 Entangled States 

The Einstein-Podolski-Rosen (EPR) [012999, 2] paradox, concerning the spatial correlations 
of two entangled quantum systems, is perhaps the most famous example of quantum dy- 
namics that is incompatible with a classical view. An entangled state is one that cannot be 
factored into the product of the individual particle wavef unctions. As a result, the state of 
one particle is necessarily correlated with the state of the other, and these correlations differ 
from those allowed by classical mechanics. Entanglement in quantum mechanics is normally 
raised to explore aspects of non-local effects and hidden variable theories. Due to the close 
proximity of nuclear spins and the fact that the ensemble is in a highly mixed state, the 
NMR measurements discussed below do not address these issues. Nevertheless, we can use 
the ability of liquid state NMR to simulate strong measurement to show that the behavior 
of an entangled state is inconsistent with a simple classical picture. 

The entangled state = -^(|00) + |H)), otherwise known as a Bell state, is given by 
the density matrix 

PBell = \ + 2I Z Sz + 2IxS X — 2IySy^j . (13) 

The above state can be prepared directly from the pseudo-pure ground state 1 00) by the 
transformation 

U = e~ iIxS ^ (14) 
which is implemented by the pulse sequence 
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(15) 
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Readout pulses can then be used to verify the creation of this Bell state, as shown in Fig |j. 

One of the advantages of working with an ensemble is that we can introduce a pseudo- 
random phase variation accross the sample to simulate the decoherence that accompanies 
strong measurement. A pseudo-random phase variation in a given basis can be achieved by 
rotating the preferred axis to the z-axis and then applying a magnetic field gradient followed 
by the inverse rotation. This leads to the pulse sequence 



[grad(z)} 



7T 



[grad(z)} 



(16) 



It can be shown that such a measurement also "collapses" the S spin along this direction. 
Thus, half the magnetization is along the +x-axis and the other half is along the -x-axis 
leaving zero magnetization in the y-z plane. This is verified in our experiment by applying 
a series of readout pulses to confirm the creation of the 2I X S X state which corresponds to 
"collapsing" the pseudo-pure Bell state along the x-axis. The experimental results are shown 
in Fig H 

An incoherent mixture of entangled states is easily generated by the pulse sequence 
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applied to p eq (Eq. yielding the reduced density matrix 
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Suppose one wishes to measure the polarization of spin I along the x-axis and spin S 
along the z-axis. One possibility is to use selective RF pulses to rotate the desired axis (x in 
this case) to the z-axis, apply a z-gradient, and then rotate back to the x-y plane to observe 
the induction signal as in Eq. [TR Alternatively, one could rotate the desired measurement 
axis of one of the spins to the z-axis, rotate the other spin to the x-y plane and then spin- 
lock the sample on resonance. In this latter case the inhomogeneities in the RF pulse and 
background field serve to effectively remove any signal perpendicular to the desired axis, and 
the induction signal is the same as in the first case. Thus for example, if a measurement 
along y for spin / and along x for spin S were required, observing the induction signal after 



the sequence 



7T 



— [spinlock]^ 



(19) 



Because one of the spins remains along the z-axis while the receiver is in phase with the 
other, the measured signals are anti-phase. The spectrographic traces shown in Figs. 6a- 
d indicate the results of the measurements Tr (4I x S y pf), Tr (4IyS x pf), Tr (4I y S y pf), and 
Tr (AI x S x pf), respectively. The traces show the Fourier-transformed induction signal read on 
the 13 C channel, with absorptive peaks in phase along either the +x- or +y— axis, depending 
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on which axis the carbon nucleus was spin-locked. Notice that Fig. ^|(d) shows the same 
anti-phase signal as the other spectra, but "flipped" by 180°. 

The results of the four plots, taken together, show a simple inconsistency compared to a 
model of only two uncorrelated classical magnetic dipoles. The product of the four traces 
has an overall factor of —1, yet each magnetic moment is measured twice so that their signals 
should cancel. Each measurement is assumed to record either the x or y polarization if each 
dipole is measured independently of the state of the other. 



15 



6 Quantum Logic Gates 

NMR provides a means whereby it is possible to analyze experiments as building blocks 
for a quantum information processor (QIP). Because spin \ particles can have two possible 
orientations (up or down), it is natural to associate spin states with computational bits. 
Further, NMR experiments can be viewed as performing computations on these quantum 
bits (qubits). 

6.1 Pulse Sequences As Logic Gates 

Suppose we wanted to implement the controlled-NOT (c-NOT, or also XOR) gate, common 
in computer science, using NMR techniques. A c-NOT gate performs a NOT operation on 
one bit, conditional on the other bit being set to 1. The action of a c-NOT gate is summarized 
by the truth table 



-A-input 


B input 
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where the True and False values have been associated with up spins and down spins, respec- 
tively. The above truth table corresponds to a unitary transformation that implements 



|00) - 


- |oo> 


|01) - 
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|10> - 


- Ill) 
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(20) 



In a weakly coupled two-spin system, a single transition can be excited via application 
of the propagator, 
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which for a perfect out = it rotation becomes (to within a phase factor) 
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(22) 



It is clear that exciting a single transition in an NMR experiment is the same as a c-NOT 
operation from computer logic. In NMR terms, the action of the c-NOT gate is to rotate one 
spin, conditional on the other spin being down. Figure [7] shows the result of performing a 
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c-NOT on p eq . While NMR is certainly capable of implementing the c-NOT operation as is 
done on a classical computer, that alone does not demonstrate any of the quantum dynamics. 
Gates implemented on a quantum information processor which have no classical counterpart 
are of much more interest. An example of such a gate is the single-spin Hadamard transform, 



i ii-.+h 



(23) 



which takes a spin from the state |0) into the state 4|(|0) + |1)). This is just a n rotation 
around the vector 45° between the x and z axes. A spectrum demonstrating the application 
of the Hadamard transform to the equilibrium state p eq is shown in figure [8]. The c-NOT 
and single-spin rotations can be combined to generate any desired unitary transformation, 
and for this reason they are referred to as a universal set of gates. [69991] 

Analysis of conventional NMR experiments in terms of quantum information processing 
has led to a great deal of insight into areas such as the dynamics of pulse sequences for logic 
gates [29992], and the effective Hamiltonian for exciting a single transition [39992]. 

6.2 The Quantum Fourier Transform 

One of the most important transformations in quantum computing is the Quantum Fourier 
Transform (QFT). The QFT is a necessary component of Shor's algorithm, which allows the 
factorization of numbers in polynomial time[49992], a task which no classical computer can 
achieve (so far as is known). Essentially, the QFT is the discrete Fourier transform which, 
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for q dimensions, is defined as follows 



QFT q \a) — > —p. exp(2niac/ q)\c) 
c =o 



(24) 



This transform measures the input amplitudes of \a) in the |c) basis. Notice how the quantum 
Fourier transform on |0) will create an equal superposition in the |c) basis, allowing for 
parallel computation. In matrix form the two-qubit QFT transformation QFT 2 , is expressed 
as 



/ \ 
1111 
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(25) 
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As formulated by Coppersmith [59992], the QFT can be constructed from two basic unitary 
operations; the Hadamard gate Hj (Eq. ^), operating on the jth qubit and the conditional 
phase transformation Bjk, acting on the jth and fcth qubits, which is given by 



Bjk 



( \ 
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10 
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^000 e^ k j 



^1(1-2^)1(1-2^) 



(26) 
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where 6jk = ^j- The two-qubit QFT, in particular, can be constructed as 



QFT 2 = F 5oi#i (27) 

The Bjk transformation can be implemented by performing the chemical shift and coupling 
transformations shown in Eq. Figure |^ shows the implementation of the QFT on a two 
spin system. The spectra show the 90° phase shifts created after the QFT application. 
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7 Conclusion 



Several basic but important concepts relevant to QIP are illustrated by experiments on a 
liquid-state ensemble NMR quantum information processor. While pure quantum mechan- 
ical states are not achievable here, the creation and application of pseudo-pure states is 
demonstrated. Tests of spinor behavior and entanglement are also described, illustrating 
quantum mechanical dynamics. Finally, building blocks (the Hadamard, c-NOT, and QFT) 
for a more complicated quantum computer are also introduced. 
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Below are the Captions for the Figures. 

1. At room temperature, the equilibrium state of chloroform molecules in solution is 
described by jjl z + 7sS , z . In order to create a pseudo-pure state from the equilibrium 
state, it is convenient to start with an equalized magnetization for I and S. Since the 
ratio of 7/ to 75 is a factor of four, then the spectra of I and S following a f pulse 
should reflect the 4:1 ratio in the peak heights (figure (a)). In order to compensate for 
the different electronics in the two channels, the gains of the channels were manually 
calibrated to produce the desired 4:1 ratio in signal intensity. After this, the pulse 
sequence discussed in the text was applied. A subsequent f read out pulse results in 
the spectrum of figure (b). The peaks have equal intensity, confirming the creation of 
the state ^^(I z + S z ). 

2. Once the pseudo-pure state p pp — I z + S z + 2I Z S Z has been prepared, we use readout 
pulses to generate a series of spectra confirming that the desired state has been created. 
This is done by applying ||^, ||^, and ||y' 5 read pulses on the pseudo-pure state. The 
results are shown in figures (a)-(c), respectively, on both the carbon and hydrogen 
channels. The signature of the appropriate terms in p pp is seen from the three sets of 
spectra generated. 

3. The state 2I Z S X correlates the spinor behavior of spin I to the reference spin S. The 
propagator U = e ~' l ^ Iy ^~ Sz ^ then rotates all the I spins coupled to the down S spins 
by the angle <p about the y-axis. Applying U to the density matrix 2I Z S X creates the 
state 2cos(<fi/2)I z S z + 2sm(<j)/2)I x S x , where only the first (antiphase) state is made 
observable by evolution under the internal Hamiltonian. When = 0, the state is of 
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course 2I Z S X , as shown in figure (a). When = 2n, this state is inverted, contrary to 
common intuition. The resulting spectrum is shown in figure (b). Only when = 4ir 
does the antiphase state return to its original state as seen in the spectrum (c). These 
spectra clearly demonstrate the spinor behaviour of spin \. 

4. The pseudo-pure Bell state, pseii = 2I Z S Z + 2I X S X — 2I y S y , created by the application 
of the propagator, U = e lIxSy7T on the pseudo-pure state discussed above can be verified 
by applying a series of readout pulses on pseii- Using the read pulses ||^, ||f , ||^, and 
|| ^ on p Beii-, figures (a)-(d), respectively, and observing the resulting spectra on both 
the I and S channels confirms both the signature and the individual terms of pBeii- 

5. We simulate a strong measurement (one that collapses the wave function along a pre- 
ferred basis or axis) on a Bell State using magnetic field gradients. An x-measurement 
on the I-spin is imitated by applying a selective x-gradient to it. Since the two spin state 
is entangled, this measurement necessarily collapses the S-spin along the x-direction. 
Subsequent measurements confirm that both I and S spins have transformed identically 
and that they are aligned along the x-axis. This was verified by observing the creation 
of the 2I X S X state where in (a) we observe immediately after the "strong measurement" 
in both channels and see zero signal as expected. In (b) we show that | pulses along 
the x-axis has no effect and in (c) we verify that a |]^ pulse indeed creates an anti- 
phase signal on the carbon channel and a pulse creates an antiphase signal on the 
hydrogen channel. 

6. Strong Measurements After EPR Preparation. All four measurements are made on 
the Carbon (1st spin) channel, and show the expected anti-phase correlation, (a) 
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Measurement of correlation I y S x , in phase along the +y direction, (b) Measurement 
of correlation I x S y , in phase along the +x direction, (c) Measurement of correlation 
IySy, in phase along the +y direction, (cf) Measurement of correlation I X S X , in phase 
along the +x direction. Note that the last spectrum is "flipped," or inverted, with 
respect to the other three. 

7. The above spectra show the implementation of a controlled-NOT (c-NOT) gate on the 
equilibrium state of 13 C-chloroform. The spectrum on the right represents the readout 
on the I spins, and the spectrum on the left is the readout on the S spins. Both spectra 
have the expected appearance and confirm the creation of the state I z + 2I Z S Z , the 
expected state after application of the c-NOT. 

8. The Hadamard gate H is a one bit gate that can be geometrically interpreted as 
a 7r rotation about the ^{x + z) axis. If the net magnetization is along the +y 
direction then the Hadamard gate should simply rotate it to the — y direction (figure 

(a) ). However, since any it rotation about an axis in the x-z plane performs the same 
transformation, H was also applied to an initial +x magnetization. The result (figure 

(b) ) shows how the magnetization was sent to the z-axis, as expected. 

9. The two-qubit QFT was implemented by applying a Hadamard gate on the first spin, 
a conditional phase operator, and a Hadamard on the second spin. A Hadamard gate 
can be performed by a simple combination of three pulses: j\ x — 7r| y — f |- x - Because 
it was performed on the thermal state, the initial Hadamard was simplified to a ir/2 y 
pulse. The conditional phase change operator, B m , was implemented by delay 2 and 
pulses 3 to 6, where pulses 4 to 6 are a 7r/4 z-rotation. The final Hadamard gate 
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was implemented by the three pulses labeled 7. The phase difference of each peak on 
the spectra shows the two-bit QFT's ability to separate input states by 90 degrees. 
After the application of the QFT, the spins were phase shifted by 45 degrees and were 
allowed to evolve for a time 1/4 J in order to bring out the phase differences. 
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Equalizing Magnetization 
Carbon Channel Hydrogen Channel 
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Equalizing Polarizations 
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Tomography of Pseudo-pure States 

Carbon Channel Hydrogen Channel 
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Creating a Pseudo-Pure State 
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Demonstration of Spinor Behavior 
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Demonstrating Spinor Behavior 
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Tomography of Bell States 
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Strong measurements on a Bell state 
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Demonstration of EPR Correlations 
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EPR Correlations 
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Demonstration of the c-NOT Gate 
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Demonstration of the Hadamard Transform 
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Demonstration of the QFT 
Carbon Channel Hydrogen Channel 




Two-spin QFT 
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